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Abstract 

We trace here instantons through the analysis of pure Yang-Mills gluon 
Green functions in the Landau gauge for a window of IR momenta (0.4 GeV 
< k < 0.9 GeV). We present lattice results that can be fitted only after 
substituting the BPST profile in the Instanton liquid model (ILM) by one 
based on the Diakonov and Petrov variational methods. This also leads us 
to gain information on the parameters of ILM. 

1 Introduction 

An appealing approach to analytically understand some of the non-perturbative 
features of QCD is the evaluation of quantum fluctuations around topologically 
non-trivial classical solutions through the expansion of the path integral around 
these solutions. In fact these considerations are often generalized to configura- 
tions which are not exact solutions of the field equations but close to them and 
which we will name quasi-classical field configurations. 

Famous examples of non-trivial solutions of classical equations of motion are 
instantons PU E] . Quasi-classical solutions considered in instanton liquid mod- 
els E] provide a successful connection between the instanton zero modes and 
the QCD chiral symmetry breaking (See ref. [5] for a good review on the subject). 

1 Unite Mixtc dc Recherche du CNRS - UMR 8627 
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More recently, it has been proven in ref. jB] that instanton model predictions 
for quark-quark interaction agree with non-perturbative QCD lattice results bet- 
ter than those from Schwinger- Dyson (SD) models with a perturbative struc- 
ture of the QCD interaction (vector quark-gluon coupling parametrization) . In 
general, a dominance of instanton-induced effects on the dynamics of the QCD 
light-quark sector seems to emerge, although it is not excluded that SD models 
with pseudo-scalar and scalar quark-gluon couplings might capture this emerging 
instanton physics. On the other hand, despite instantons have been first consid- 
ered as a possible explanation of confinement [Zj, it is now generally accepted 
that they do not generate the area law for Wilson loops 2 . 

We have recently argued [S] that instantons, or instanton-like structures, have 
dramatic effects on the low momentum Green function in Yang-Mills theories and 
that they can explain the observed oc k A behaviour of the non-perturbative MOM 
QCD coupling constant computed on the lattice. Notice that this remark, not 
only advocates in favor of the presence of these quasi-classical structures in the 
lattice gauge configurations, but also indicates that the quantum fluctuations do 
not contribute significantly to the Green functions in this momentum regime. In 
the present paper we try a low momentum description of two- and three-gluon 
Green functions through Instanton liquid model (ILM). 

The succesful description of the MOM QCD coupling constant in the low- 
momentum regime is based on the sum-ansatz approach, that builds the classical 
solution as a linear combination of modified instantons. Although instanton 
profile modifications play no role in obtaining the coupling constant, at least in 
the first approximation, any description of two- and three-gluon Green functions 
makes mandatory to further elaborate on the nature of this modifications. To 
this goal we will follow the Diakonov & Petrov (DP) sum-ansatz approach pQ. As 
will be discussed later, different aspects of this approach have been criticised by 
Shuryak and Verbaarshot, but it provides us with a framework able to estimate 
instanton effects through analytical or semi-analytical computations which seem 
to work reasonably well. This is the "phenomenological" point of view which 
we adopt in this paper to extract some understanding of the low-momentum 
behaviour of lattice Green functions. 

The paper is organized in six sections. In section 2 we discuss the pattern of 
the running with momenta of the QCD coupling constant. Section 3 is devoted 
to study the instanton profile modification within the DP approach. In section 
4, we show that the low momentum behaviour of lattice gluon Green functions is 
rather well described by ILM only after including instanton profile modifications 
and discuss how the large instanton density obtained from the fits is expected 
to be reduced by light dynamical quarks in full QCD. In section 5, the effect of 
instanton radius distribution is discussed. We finally conclude in section 6. 

2 Another set of solutions of classical field equations, merons UJ, has been recently reexamined 
as candidates to explain confinement |10j . 
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At the end of the day, of course, rather important questions still remain open: 
why do quantum effects appear to be suppressed in this low momentum range ? 
Has this some connection with confinement ? 

2 The QCD coupling constant: four regimes 

We have presented in Ref. jS] a preliminary claim of instanton dominance at 
low energy by analyzing in Landau gauge the following ratio of pure Yang-Mills 
Green functions: 

4vr (G( 2 )(F)) 3 ' 

which is a non-perturbative MOM definition of the coupling constant, where 
is the gluon n-point correlation function and Z M q M = k 2 G^ is the gluon 
propagator renormalization constant in MOM scheme 3 . 

In Fig. ^we show on a log-log plot, that a roughly-/c 4 power law is satisfied 
by the lattice evaluations of a s (k 2 ), eq. ((H), up to 0.8-0.9 GeV, for three different 
lattices, strongly supporting a quasi-classical description [Hj. 

A very striking feature of the results shown in Fig. Q obtained with a low 
statistics of 20 configurations, is that there is a rather sharp transition at ~ lGeV 
between two regimes: below this scale a s (k 2 ) does not seem to fluctuate much, it 
agrees well, with a small \ 2 i with the expected k A linear behaviour in spite of the 
small statistics and it complies with an instanton-like picture, while above that 
scale the data suddenly deviate from the k A law and become apparently fuzzier 4 , 
a fuzziness which is confirmed by the larger x 2 which affects any smooth fit. A 
tempting interpretation is that this fuzziness has to do with a strong influence of 
quantum fluctuations which increase suddenly above 1 GeV. We may understand 
this as follows: for a given profile and a given instanton radius, the ILM predicts 
no statistical fluctuation of the Green functions. The average over the instantons 
locations and their color orientation does not create any noise on Green functions 
which are translational and color rotation invariant. Only the dispersion (to be 
studied at length in this paper) of the instanton radius as well as varying effects of 
the neighbouring instantons produces some statistical noise. On the contrary the 
quantum contributions to Green function are generated by statistical fluctuations 
of the gauge fields around zero and Green functions appear as correlations in this 
statistical system. Some confirmation of this interpretation results of the analysis 
of the same gluon Green functions after applying a cooling procedure that kills 
short- distance (quantum) correlations [TT] . 

Altogether we may distinguish four regimes: 

3 Wc keep the name "coupling constant" for this well defined quantity although it could be 
argued that this name is not really appropriate in the low momentum regime. 

4 The smooth curve shown in fig 1(b) in [5j has been reached with 1000 configurations. 
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• Above 2.6 GeV we have shown ^21 EH] that the lattice data were domi- 
nated by perturbative QCD with a significant non-perturbative correction 
describable via OPE by the expectation value of A 2 . This means clearly a 
dominance of quantum fluctuations with small corrections from the A 2 con- 
densate which may be generated by the quasi-classical solutions. Indeed the 
quasi-classical solutions, being large structures, are seen by the hard propa- 
gating gluons as an effectively translational invariant background. It is then 
easy to show that their dominant effect is amenable to an OPE treatment 
of the lowest dimension operator: A 2 |14j . 

• Between ~ 0.4 and ~ 0.9 GeV the quasi-classical contributions dominate 
and the quantum effects are strongly depressed. 

• The 1.0 to 2.6 GeV region shows a strong quantum effect. However, it is 
not at all describable in terms of perturbation theory. A description in 
terms of quantum fluctuations in a quasi-classical background should be 
tried. The latter background can no longer be treated as simply as in the 
large momentum regime. Other non-perturbative effects may also play a 
role, for example related to confinement. 

• The very low momentum region below ~ 0.4 GeV is still "terra incog- 
nita" and, being of the order of 2n/L, (L being the lattice length) possibly 
strongly sensitive to finite volume artifacts. 



3 A modified profile for the instanton-liquid model 

In the following, we will try to describe the Green functions assuming that the rel- 
evant quasi-classical solutions are instanton liquids. Let us compute the quantity 
in eq. (JTJ). We describe gluon fields in the low energy regime as a superposition 
of modified instantons, a quasi-classical background dominance being assumed. 
In this framework, the gauge field in Landau gauge will be given by: 

J B;(x) = 2^fl S ^^^^(J^l) , (2) 

where z % (pj) are the center (radius) of the instantons, 77" is known as 't Hooft 
symbol, Ry? are color rotations embedding the canonical 577(2) instanton into 
the 577(3) gauge group, a = 1, 3 (a = 1, 8) is an SU(2) (577(3)) color index, and 
the sum is extended over instantons and anti-instantons. 

The classical solution for an isolated instanton is the standard BPST one, 
0(£) = l/(£ 2 + 1) P- Nevertheless, the superposition (j2J) with a BPST profile 
is not a solution of Yang-Mills equations (since they are not linear). Hence, 
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Figure 1: The low momentum k behavior of a s (k 2 ) on log-log plots for three different 
lattices ((5.6,24), (5-7,32), (5.8,24)) using 20 configurations. Notice that above ~ 1 
Ge V, no smooth line can be easily drawn joining the lattice data. 



we will discuss below a parameterization of the profile function inspired by an 
approximated minimization of the action for a finite density of instantons 5 . 

Then, just by assuming random color orientation and instanton position 6 , 
from eq. (J2|) we obtain (see [£]): 



g^ ) {e)= l -< P H(kpf> 
g^ ) {k\k\k*) = ^- k <p*I{kpf> (3) 

where n stands for the instanton density, < ■ • • > denotes average over the 
instanton radius for a given radius distribution p(p), normalized to 1, and where 



I(s) = / zdzJ 2 {sz)(j)(z) , s > ; (4) 

s Jo 

J2 being the second order Bessel J function. The factor g n for n-points Green 
functions in l.h.s. of eq. (J3J) comes from the factor g in the l.h.s. of eq. (J2J). 

5 The function cj> in eq. (J5J) is related to the function P in eq. (6) of ref. |B] by the relation 

p(e) = zm/e- 

6 This means that we neglect the color correlation which might exist, for example, between 
neighbouring instantons, an assumption which is usually done and which amounts to consider 
this instanton liquid as not being ordered. 
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Then, from Eq. JI} we get: 

as ^ k) = 18^ < pH{k P f >3 ' (5) 

As a first approximation, if we consider all instantons to have the same radius, 
we obviously obtain: 

< P 9 Hkp) 3 > 2 = , . 

< p 6 /(A;p) 2 > 3 ' 1 ' 

and recover tin GXcict k - power law for any instanton profile. 

In © the influence of the profile will only appear as a sub-leading contribu- 
tion, that will also depend on the instanton radius distribution (See sectional), 
while in (J3J) the leading contributions to the Green functions depend on the profile 
and will therefore be used by us to gain some understanding about the instanton 
radial shape and radii. Indeed, it will be manifest (see fig. EJ) that the BPST 
profile cannot account for the low momentum behaviour of Green functions. 

The variational Diakonov & Petrov equation. 

If the QCD vacuum can be understood as an instanton liquid of finite density 
(See jHj, for example), the BPST profile will no longer be valid (it is only a zero 
density limit) especially at large distance from the instanton center where the 
overlap of neighboring instantons becomes important. Being interested in the 
low momentum regime we cannot neglect these effects. 

A possible method to include the effect of instanton interactions, is to study 
the profile that minimizes the action of the instanton ensemble. Such a procedure 
was used in [3], where, through the Feynman variational principle, the equation 



(7) 

was obtained for the best profile, 0, eq. (6.7) in jlj 7 . In the last term (3{p) = 
87i 2 /g 2 (p) where g(p) is the running coupling constant evaluated at the instanton 
radius scale and Cn c is a factor containing the quantum corrections that multiply 
the instanton liquid partition function and are basically defined by the functional 
determinants in eq. (2.19) of jlJ.The value of «dp which represents an average 
classical effect of the other instantons on one of them is given by 

a l F = l 7o 2 n§4 < /"°°^ 2 2 (-)> ■ (8) 



3' U M Jo \P 



7 We correct for a factor x 2 in the last term which misprinted in the quoted paper. The same 
for the square of the function <f> in eq. |JSJ. 
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where the parameter 70 is related to the instanton interaction strength, 



7 2 - 27A ^ 2 (9) 

Thus, if we assume the same radius for all the instantons, this equation reduces 
to 

. 2 ^ 9n c k\ r,w* 



a ^w^ n k dx * {-p) • (10) 

The last term in eq. (J7J), involving a functional derivative of the functional 
determinants, is hard to compute and explicitely violates scale invariance (in the 
sense that </>(-) is still a solution if n is divided by a 4 ). Let us remark that this 
is true only if we do not take into account the p dependance of «dp due to /3(p). 

Now assuming that «dp is approximately constant, it is argued in jlj (dis- 
cussion before eq. (6.1)) that one can neglect this term to get a scale invariant 
equation governing the large distance shape of the profile function assumed to 
be dominated by classical interaction. Of course, we should borrow from the ne- 
glected term some scale invariance breaking to fix the instanton size and match 
to the BPST solution at small distances. This can be done e.g. by explicitely 
writing the profile as a function of \x\/p, p defining the instanton size through 
the condition 0(1) = 1/2. 

The term in ojdp enforces a squeezing of the instantons and therefore controls 
the large x behaviour of cj)(x). For \x\ > p we can neglect in Eq. (J7|) the non 
linear terms obtaining the following Bessel equation: 

^ + l 1 ' l dR-( 4 + a - 1 ' 2 )} *{t) = °' (11 > 

Therefore we have: 



U ~ cK 2 (a DF \x\) ~ c^=^, (12) 
P J x > p vfI 

where c is an unknown coefficient which could be determined from the scale 
invariance breaking condition. Let us make further remarks about Eq. (J2J) after 
neglecting the Cn c term: the authors of [3] claim that at small x the term 
proportional to a^pX 2 becomes small and therefore can be neglected, recovering 
the instanton equation. This term is nevertheless not negligeable wih respect to 
cf) — 3<p 2 + 20 3 which tends to when x — > 0. Furthermore this equation has no 
solution going to at infinity and to 1 at x — because a singularity emerges at 
some small finite value of x. 
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For the above reasons, we prefere to use in the next sections instead of direct 
solutions of eq. the following parametrization: 



x\\ _ Qppp) 2 K 2 (appx) . . 

P J 2 t , (a DP pf v , ' [6) 

1 H K 2 {a D px) 

that behaves as ^(ccdpI^I) at x ^> p 8 and at short distances as: 



^ — =^—^ + 0(a DP x) , (14) 
VP/ P + x 

where for arjpx C 1 we indeed recover the BPST instanton. Remember that 
the zero «dp limit is also, from eq. (|10p. the zero density limit. Notice that 
the imposed constraint, expressed by eq. (JHJ), of a behaviour a la BPST at 
«dp x <C 1 is de facto the scale invariance breaking condition in our approach. It 
differs slightly from DP's proposal 0(1) = 1/2 since 0(1) = 1/2 — («dpp) 2 /16 + 
0((q;dpp)4). For the values of odp and p which we will use to fit the lattice 
data both are in practice equivalent. We show in fig. |2]a plot of this profile for 
x/p > 1/2 compared to a numerical solution of Eq. (JJJ) equal to 1/2 at x/p — 1 
and going to at infinity (for «dpP = 0.675). In the following, we will consider 
eq. (|13|) as our optimal choice for the profile function. 



4 Green functions 

The goal of this section is to apply the parameterization in eq. ([13)1 to fit our 
numerical results, obtaining thus p, q?dp and the instanton density. We start by 
writing the gauge field as the addition of a classical part, -B^(x), and a quantum 
one, Qb(x), depending in general on the classical background: 

^(x) = B«(x) + (Q B );(x) , (15) 

then, for the two-point gluon Green function, we can write 9 : 

(^(Q)4(x)) = (ggo )^(x)) + ((g B )*(o)(g fl )t(x)> , (i6) 

In ref. |15| . in the context of a constrained instanton model, the authors propose Ansatze 
to account for large-scale vacuum field fluctuations also by matching similarly large and short 
distances limits of their constrained instanton equation. They also obtain solutions decaying 
exponentially at large distances. 

9 Crossed terms vanish if the background corresponds to a local minimum of the action i.e. 
to a classical solution. 
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Figure 2: Comparison between our parametrization of the profile function wit arjpp = 
0.68 and the exact solution of the Diakonov and Petrov equation. 



or after Fourier transformation, 

c2t ice (fc 2 ) = Gg(fc 2 ) + G%\k 2 ). (17) 

The working hypothesis we derive from the interpretation of fig. ^ is that in 
a certain region of momenta, say below ~ 1 GeV, quantum effects are strongly 
suppressed and only classical properties are seen. Of course, we introduce this 
hypothesis based on a phenomenological observation, but it is not unconceivable 
that an intrinsically non-perturbative phenomenon like confinement could cause 
the disappearance of quantum correlations at distances larger than the confining 
scale. 

4.1 Estimated corrections to the instanton liquid model 

From our data we have been led to fit the bare Green functions to the r.h.s of 
eq. (J3J) in a range k m i n < k < k mSuX without any correction. It is far from obvious 
that we can neglect either quantum corrections or lattice truncations of the quasi- 
classical model. This subsection is devoted to justify our choice. For the sake of 
clarity, we will pay the price of anticipating on the results of some fits which will 
be later detailed. 

(i) The Green functions G^ n > estimated from three different lattice spacings 
match reasonably to each other, after a mere rescaling to a common value for 
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Figure 3: (a) We plot the gluon propagator, k 2 G^ 2 \ for the three lattices after a 
rescaling that matches data from (3 = 5.7 and (3 = 5.8 to those from (3 = 5.6. The dotted 
line is obtained from the BPST profile with p = 0.3 fm and the solid line is the best fit 
using our profile parameterization for the gluon propagator. BPST cannot manifestly 
describe lattice data, (b) The points are lattice evaluations of symmetric vertex, k^G 1 ^' . 
The solid line is the best fit for the vertex with our profile parameterization, the density 
being required to be the same for both propagator and vertex, and the dotted is the best 
fit with a free density in the fitting. 



some fc max of the order of 1 GeV (see Fig. EJ). These matching coefficients are 
close to 1, approximately in the ratios 0.95 / 1.0/ 1.05 for (3 = 5.6,5.7,5.8. 

(ii) It is difficult to know precisely which mechanism drives these matching 
coefficients slightly away from 1: quantum corrections, ultraviolet /infrared 
cutoffs on quasi-classical solutions. It is also impossible to know if these cor- 
rections are multiplicative, additive or of a more complex nature. The main 
lesson is that these corrections are small, and we decide for convenience to 
describe them by a multiplicative rescaling factor. 

(iii) This factor is defined as 

M m ) (U2 n -i\ 

rw^,,,- 1 ) = latti "T — , (is) 

G$(k 2 ) 

where a -1 is the regularization scale i.e. the inverse lattice spacing and 
G$(k 2 ),m = 2, 3 are given in eq. ©. The functions T^(k 2 , a" 1 ) n/3/6 10 , 
n being the instanton density, are plotted up to an unknown global constant 
in fig. Ufa). G^(k 2 ) is here the best fit to be discussed later. These plots 
show the good matching of different lattice spacings after performing a 



3 The factor (3/6 = 1/g 2 comes from the factor g 2 in eq. (J3J). 
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constant multiplicative rescaling (figHJb) shows the result of this rescaling). 
They show a wide parabola having its minimum around k = 0.65 GeV, this 
minimum is around 7 to 8 which corresponds to the order of magnitude of 
the instanton density we will derive (the other factors being close to 1). 

One possible explanation of this parabolic behaviour is an expected increase 
of quantum corrections towards larger k and, towards lower k, a relative 
increase of the ratio due to the fast decrease of the denominator in eq. (fTHjl : 
for example the fast increase at low k might be due to additive quantum 
corrections which are only visible when the quasi-classical background is 
very small. Additive corrections are anyhow necessary at k = since the 
lattice data are non vanishing while the denominator of eq. (|18|) is zero. 
Whether this non vanishing of the lattice propagator at k = is a finite 
volume effect violating Zwanziger's theorem [TZj, will not be discussed in 
this paper since, as already mentioned, we restrain from discussing the 
finite- volume sensitive region below ~ 0.4 GeV. 




0,4 0,6 0,8 1 0,4 0,5 0,6 0,7 0,8 0,9 

(a) (b) 

Figure 4: (a) We plot the ratios nT ( - 2 \k 2 ,a~ 2 )P/Q where n is the instanton density 
and I^ 2 ) defined in eq. U8\) . for the three lattices, the squares, triangles and circles 
corresponding respectiveley to (3 = 5.8, [3 = 5.7 and (3 = 5.6. The horizontal axis 
represents the momentum in GeV. (b) The same points as in (a) but with the mul- 
tiplicative rescaling performed. Different (3' s coincide quite well, all show a parabolic 
type dependence in k. 

(iv) If we stick to the momentum range plotted in fig. IHb) these ratios do not 
deviate from a constant by more than a few percent. We will therefore per- 
form our fits only in this range k m i n — k max 11 and approximate T^(k 2 , a -1 ) 

11 To be precise our fits use the window fc min = 0.44 and /c max = 0.89. 



11 



by a constant: 



T^(k 2 ,a^) ~ T^ia- 1 ) , (19) 

L.2 < t,2 < u2 
min~ n ' ~""max 

(v) Since from fig. 0Ja) the different r^ 2 ^(a _1 ) differ also by no more than a 
few percent we decide to take from now on all these r^(a _1 )'s as equal at 
the cost of an expected discrepancy between the fitted densities of a few 
percent. We take 

Va"\ I^V 1 ) = r< 2 ) {a~\5.6)) . (20) 

where a -1 (5.6) is the inverse lattice spacing for f3 = 5.6. 

(vi) We now wonder if the hypothesis eq. (|2U|) can be extended to the coefficients 
T^ 3 \a~ 1 ). The small dependence on k and on a" 1 can be seen, although 
with larger errors, in fig. Efb) where the different lattice data seem to fit 
one common solid curve representing the model. The next question is the 
validity of the hypothesis 

r^a- 1 ) =r( 2 V 1 ). (21) 

To test the hypothesis I^ 3 ) ~ we can look at fig. El The vertical 
axis corresponds to the ratio (a _1 (5.6)/r^ 2 ^(a _1 (5.6)) or equivalently 
to 7i® /rS 2 ' where we define rv~ m ' = nr( m )(a _1 (5.6)), n being the instanton 
density. The horizontal axis corresponds to the instanton radius p, supposed 
to be the same for all instantons. For that instanton radius a best fit of the 
bare Green functions is performed according to the formula: 

Gl:L c (k\a~i) = G^(k 2 ) r^(a- 1 (5.6)) (22) 

with G^(k 2 ) given by eq. 0- The plot shows for each p an agree- 
ment between different /5's which is a surprise: it tells that the ratio 
r( 3 )(a _1 )/T( 2 )(a -1 ) is almost independent of the lattice spacing even though 
it does strongly depend on p. The dotted curve shows the \ 2 of the 
common fit to and G^ 3 \ The smallest x 2 corresponds to a ratio 
(a _1 )/r^ 2 - ) (a~ 1 ) ranging from 0.8 to 0.9, i.e. close to 1. This also cor- 
responds to a value of the radius around 0.3 fm which is the one favored by 
phenomenology. This leads us to consider eq. (|2*Tf as quite reasonable. 

(vii) All the arguments up to now have shown an approximate equality to a few 
percent of all r( m )(a" 1 )'s. We are left with an unknown global constant. 
The fact that these factors are so close to each other suggests that the 
corrections to the instantonic contribution (whatever their origin may be) 
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are small. It is then natural to expect that these correction factors are also 
close to 1. This is confirmed by a result which will be detailed at the end of 
this paper via the following argument: if we initially take T^ 2 ' 3 ^ ^ 1, apply 
MOM renormalisation and then eq. ([I]), we obtain: 



«Latt(^) = 




where the term above the curly bracket in the r.h.s. introduces a supposed- 
to-be-sub-leading correction, depending on the regularization parameter, 
that we neglect in our analysis for the present work (and did the same 
in jH]). Let us anticipate some of our next results: the reasonable agreement 
between the instanton density obtained from the coupling constant ctLatt(fc) 
in table 01 with the one in table Q implies a ratio (r^) 2 /(r^) 3 close to 
1. The latter, combined with /T^ ~ 1 presented in the preceding 
item, ends up with an approximate justification of our hypothesis that all 
r( n )(a _1 ), n = 2,3 ~ 1 for all (3's in the range considered here. 

Thus, our final working hypothesis is: 

r^(fc 2 ,a- 1 ) = l for m = 2, 3 and k min < k < fc max . (24) 

with k min = 0.44 GeV and A; max = 0.89 GeV. 



4.2 Numerical results. 



The numerical data that we shall exploit systematicaly all over the paper result 
from two simulations on a 24 4 lattice with bare coupling constants given by 
(3 = 5.6 and (3 = 5.8, and a simulation on a 32 4 lattice for (3 = 5.7 12 , in the 
Landau gauge. The gluon correlations functions are obtained from only 20 field 
configurations, that are enough to manifest the classical background effects we 
search for. However, more field configurations are available for some of the lattices 
and we used them to check that those exploited are not particularly biased 13 . We 
calibrate all these simulations with the ratios of lattice spacings for different /3's 
given in ref. [TH] and a^attiP = 6-0) = 1.97 GeV (this last value was used in 
ref. ^1 consistently with the very precise measurement of the lattice spacing 
resulting from a non-perturbatively improved action in ref. [20j). 

Coefficients in ^H] are fitted for /3 = Q/g 2 > Q. In this work, we make use 
of rather low values (3 < 6.0, in order to have larger volumes. These low values 

12 With the two lattices for (3 = 5.6, 5.7 we simulate practically the same physical volume and 
are thus in position to check lattice spacing artifacts. 

13 Our Jacknife errors estimate also the dependence on the set of configurations. 
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Figure 5: We plot the ratios /n^ where we define = rar( m )(a -1 (5.6), n being 
the instanton density. These are fitted from two-point and three-point Green functions 
for a fixed value of instanton radius, p, represented on the horizontal axis. The dotted 
line joins the optimal x 2 /d.o.f . (on left y-axis) as a function of the radius. 



of (5 show for the quantity a s a good scaling with results at f3 > 6.0, a signal 
that these lower /3's provide reasonable results. The risk with these simulations 
is the extrapolation needed to calibrate the lattice spacing, that might be a non- 
negligible source of systematic error in our measures. 



p (fm) 


0.2 


0.25 


0.3 


0.35 


0.4 


n 

(fm" 4 ) 


P = 5.6 


26.5(1) 


13.18(5) 


7.75(3) 


5.12(2) 


3.63(1) 


P = 5.7 


27.2(1) 


13.54(5) 


7.98(3) 


5.28(3) 


3.76(2) 


P = 5.8 


27.8(2) 


13.9(1) 


8.20(7) 


5.41(4) 


3.86(3) 


a DP p 


0.393(2) 


0.527(2) 


0.675(2) 


0.836(3) 


1.005(3) 


X 2 /d.o.f. 


6.1 


3.3 


2.3 


3.3 


6.0 



Table 1: Best-fit parameters for several fixed instantons radii. The two- and three- 
point Green functions are fitted simultaneously with the some instanton density. Note 
that the stated errors are only statistical. 



We collect the results of our fits for several fixed instanton radii in tab. H An 
instanton density, although the same for both Green functions 14 , is independently 

14 Notice the difference with the fits presented in fig. [S] where different densities are assumed 
for G (2 > and G (3) . 
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fitted for each particular lattice spacing. We look thus for a remnant of a sub- 
leading dependence on the regularization parameter. The careful reader may 
have noticed that the relative density splitting between different /3's is slightly 
smaller (3 % for 5(3 = 0.1) than the splitting between the linear fits in fig. |U(4.5 
%). This is due the factor g 2 in eq. © (about 1.5 %). 

Crudeness of our quasi-classical approximation. 

It is worth to insist that our fit rely on the crude approximation given by eq. 
(124}) which is not valid beyond a few percents as discussed at length in subsection 
14.11 This is the origin of the difference between the densities in tab. H This 
is also the reason why the minimal y^/d.o.f. in table Q is of the order of two 
because we force the overall factor to be the same in fits for both two- and three- 
point Green functions. This tells about the crudeness of completely neglecting 
the sub-leading contributions from quantum fluctuations and lattice truncation 
of the classical solutions. Much better matchings might have been obtained had 
we relaxed the constraint eq. (}2"4*j) . We refrained from doing so because it would 
have needed additional input about quantum fluctuations and lattice truncation 
of classical solutions which would have been mere guesses and would not have 
yielded any stronger evidence. 

At the present stage of the work, the best we can do is to assume eq. (}2"4*|) and 
get what we believe to be a fairly coherent description of our whole set of lattice 
data. The best-fit parameterization in Tab. ^ is the best we can achieve, and 
it is not that bad, although we know that the density obtained there is affected 
roughly by a 20% of systematic uncertainty 15 . 

4.3 QCD versus pure gluon-dynamics 

Our measured densities in tab. 121 while in the ballpark of lattice estimations, 
give a contribution to the gluon condensate {g 2 G a ^ v 2 ) ~ 4GeV 4 ~ (1400 MeV) 4 
for n ~ 8 fm -4 , significantly larger than other estimates, based on QCD sum 
rules jSDj, (g 2 G a ^ 2 ) ~ 0.5 GeV 4 ~ (840MeV) 4 , or more recently |2|, (g 2 G a ^ 2 ) ~ 
0.9GeV 4 ~ (970MeV) 4 . A significant difference between the pure Yang-Mills con- 
densate and the one from QCD sum rules is however expected: in ref. [22] (see 
eq. (106) ) it is argued that light dynamical quarks might reduce the G 2 conden- 
sate by a factor 2 to 3. In fact, it can be argued that the low eigenvalues of the 
fermion determinant reduce the instanton contribution from the action and hence 
the resulting instanton density. In addition, some re-shaping of instanton profiles 
for a strongly correlated vacuum populated by fermions cannot be discarded. As 
a result, estimates in pure gluon-dynamics should differ substantially from those 

15 This uncertainty is estimated from the n^ 3 ' /n^ 2 ' ratio at the minimum x 2 m fig- HI ranging 
around 0.8-0.9 instead of 1. 
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in QCD and only the analysis of unquenched (TV/ 7^ 0) lattices simulations of 
gauge fields could help us to quantify this discrepancy. 

Not only the same analysis of this paper for unquenched simulations could 
quantify the discrepancy, but we may also appeal to the connection of A 2 - 
condensate and the ILM approach: the gluon condensate of dimension two can 
be computed by extracting the momentum power OPE correction to the per- 
turbative formula of Landau-gauge Green functions from their lattice estimates 
above ~ 2.6 GeV (see discusion in section [121 • This condensate is to be 
computed at a given renormalization point, /io, lying on the momentun range of 
OPE-dominance. 

Now, if we run the renormalization point /io of the A 2 -condensate down to 
some instanton scale, /i/, with the help of the one-loop renormalisation goup 
equations (RGE), we obtain [Tlj : 

«,) . = (A\ m ) * (Al 4 J (, + I i^^y) ; (25) 

where we match the OPE+RGE estimate to an ILM semiclassical one. Of course, 
the latter is a purely semiclassical estimate deprived of the UV-fluctuations 
around the classical minima which we assume to be subleading at the appropri- 
ated instanton scales (see again the discussion in section [2J). Thus, the instanton 
density can be roughly estimated from the dimension-two gluon-condensate. 

Although the gluon condensate for quenched simulations has been largely 
studied ^2] in literature, unfortunatly only a very preliminar analysis of the 
MOM QCD coupling constant from a lattice simulation with two dynamical-quark 
flavours is available [T^J. Furthermore, this analysis is performed for such high 
quark masses that no effect of instanton suppression is to be expected and indeed 
no significant effect has yet been seen. Progress in the unquenched determination 
of the gluon condensate will be, of course, welcome. 

Moreover, the QCD sum rules are not so accurate and neither are our den- 
sity estimates which, for instance, would be modified by the presence of other 
"instanton-like structures" or instantons deformations. 

Altogether our density estimate is close to the maximum acceptable: with a 
density of 8 per fm 4 the average distance between two neighbouring instanton 
centers is of the order of 0.6 fm i.e. twice the average radius, a really dense 
packing. 

5 The effect of the instanton radius dispersion 

We have performed satisfactory and independent fits for the three different lat- 
tices from a s to the /c 4 -formula, eqs. (0 El); and obtain the estimates for the 
instanton density shown in the second column of Tab. [21 They differ by about a 
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30 % from the profile dependent ones obtained from and presented in table 
^ The fit of a s is done on a rather firm ground since, contrarily to the fits of 
, the prediction does not at all depend on the profile neither on the radius 
but only on eq.fjBJ), of course on the hypothesis that the radius distribution is a 
delta function i.e. that all the radii are equal. a s is therefore the best quantity 
to estimate the corrections to this rather drastic hypothesis. If the k A law does 
not depend on the instanton profiles but only on the 5-function distribution, the 
corrections do 16 . 

First, we face the problem of computing the corrections to the previous 5- 
function result from a more general optics: we just consider some small width for 
the distribution, 5p, and compute in perturbations of Sp/p around the average 
radius, p. We can write: 



< p 3m I m (kp) > = I m (kp) < p 3m > + —I m (kp) < p 3m (p -p)> 

dp 

+ \j? im{m <P 3m (p-p) 2> +■■■ ■ ( 26 ) 

we take p = ~p + 5p , assume the distribution to be symmetric (at least for small 
perturbations) around the peak, i.e. < 5p >= and then obtain: 



n( n - 1) V n ( fy 4 
2 f +U {p 4 

<(p + 5 P r5p>= np n+l6 4- +o( 6 ^ 

P V P 

<(p + SpYS? > = p"+ 2 ^ + O (|Q . (27) 

Thus, we derive for the QCD MOM coupling constant defined in eq. (0) the 
following expression: 



a(k) = + Y ( 12 + 1 W(^) + 3 P 2 / 2 (^))+°(^)) > ( 28 ) 

where all the dependence on the profile function comes through the function / 
defined as 



f(kp) 



d 

dp 



In 



zdzJ 2 (kpz)(j)(z) 



o 



(29) 



16 Obviously the influence of the radius distribution fJi(p) will not be independent of the profile 
(Remember (J3J). 
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In particular, for the BPST profile, we find: 



I - (K.ikp) + K 3 (kp)) 
pf(kp) = -2 p=2 • (30) 

i - Mm 

This function ~pf(k~p) takes only negatives values for all k~p. The same is obtained 
by applying the DP-inspired profile parametrization defined in eq. (J 13)) (where 
we take ap = 0.68 from the previous section) to eq. (j2~9"|) and, as can be seen 
in plot El the functions / obtained from both profiles evolve rather close to each 
other for all k~p. This suggests a pattern that eq. (|2*H|) obeys in rather general 
manner 17 : the fc 4 -power law is corrected by a positive fcp-independent term and 
by a negative one depending on k~p\ the first one gives some pre-factor larger than 
1 to the k 4 formula and the second can be numerically simulated rather well over 
the window (0.4 — 0.9) GeV by reducing 18 the power on k, 



a(k) ~ — — (kp) 4 ~ £ . (31) 
iQimp 

This last result, eq. (|31|). altogether with the constraints c > 1 and e > 0, are 
the main and more general result of the section. The values of e and c depend 
on the particular profile and on the width of the radii distribution, as given by 
eqs. (128129)) in the small width limit. 

We will now proceed as follows. We will first estimate the parameter e by 
fitting our lattice data for MOM QCD coupling constant to the formula eq. (|31|). 
where n/c is to be taken as a free parameter to be tuned. Then, we will estimate 
the width and the instanton density, n, with the help of eqs. (|28I29|) . i.e. in the 
small width approximation. Finally, the consequences of employing more realistic 
radii distributions will be preliminarily discussed, details being in appendix iBl 



Lattice results 

In this subsection, we first collected in tab. |2]the results for the density obtained 
through the fit of the MOM QCD coupling constant defined in eq. (JIJ to the 
fc 4 -power law in eqs. ()5I6|) . for our three lattice data set. 

Then, we include in tab. |3]the results for the best fits of the same QCD coupling 
constant lattice data to the power law formula of eq. (|31|) . 

We use for both the same fitting window, 0.44 < k < 0.89, used in section 0] 
to obtain the best profile parameterization. The x 2 /d.o.f.'s in tab. El are tiny. 

17 see appendix 1X1 

18 It is in fact a matter of power reduction because k~p > 1 in the most and more numerically 
relevant part of the fitting interval, for p = 1.5 GeV -1 . 
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Figure 6: The function pf(kp) defined in the text in terms of kp for both DP-inspired 
and BPST profiles. 



Lattice 


n (fm 4 ) 


5.6(24 4 ) 


5.2(6) 


5.7(32 4 ) 


6.7(4) 


5.8(32 4 ) 


6(1) 



Table 2: Results for the density obtained from fitting lattice coupling constant to eqs. 
1516)) i.e. assuming a k 4 behavior. Errors quoted are only statistical and computed by 
the jackknife method 



This may be due to: (i) the strong correlation of our lattice estimates of the 
coupling for the different momenta, all computed from the same set of gauge 
field configurations 19 ; (ii) the small number of points in our fitting window (only 
3 for 5.8 (24 4 )). 

The best-fit parameters computed lattice by lattice in tab. H3 manifest no 
appreciable systematic deviation 20 . A global fit for the whole set seems thus to 
be appropriate. The resulting power deduced from a global fit is: 

3.91 ±0.45 , (32) 

19 Still, our jackknife estimate of the statistical error is reliable because, if the global average 
for the different momenta are correlated, the averages for different sets of gauge configurations 
are not. 

2 As expected since we found the sub-leading lattice-spacing dependent corrections to be 
small in eq. I)23|). 
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Lattice 


- (fin" 4 ) 


e 


xVd.o.f. 


5.6(24 4 ) 


6(2) 


0.2(4) 


0.068 


5.7(32 4 ) 


7(3) 


0.1(5) 


0.36 


5.8(24 4 ) 


5(1) 


0.2(6) 


0.007 



Table 3: The same that in tab. El but here fitting to eq. (jSU for p = 1.5 GeV~b The 
third row shows the x 2 /d.o.f. The fitted powers are then lower than 4 but the errors 
show that the power 4 is not excluded by our fits. 



which is compatible with both the general constraints > and also with a k 4 -power 
law. The xVd.o.f. is 0.39 for this global fit. In ref. [Hj, the same global fit was 
performed using also a few other lattices (we added estimates of a s over smaller 
volumes and larger /3's: 5.7(24 4 ), 5.9(24 4 ), 6.0(24 4 ,16 4 ); these, having only one 
point inside the fitting window, have been discarded from the present work. The 
result of the latter analysis was 3.82(8) for the fitted power. The error was there 
clearly underestimated 21 but the evaluations of the coupling constant from larger 
/5's do not practically affect the central value. This fact indicates that exploiting 
the small /3's used in the present work does not seem to introduce any sizable 
bias, at least over the low- momentum region, on the determination of a s . 



Instanton density in the small width approximation 

The global fit for the whole set of lattice data yields e ~ 0.1 for the central 
value estimate of that parameter. If we apply the exact result for / in eq. (l2~9l 
to eq. (|28|) and use the latter to estimate Sp through its best numerical matching 
to the power formula eq. (J31|) with e = 0.1, we will obtain: 



P 

This results implies that 



Sp 2 , ^ 

-V ~ 0.01 . (33) 



c~ 1 + 12 -V ~ 1.12 . (34) 
P 

Thus, densities obtained from tab. El will be about 10 % larger than those in tab. 
Elbecause of the factor c in front of the power of k. If we compute numerically the 
function / for the DP-inspired profile and follow the same procedure to estimate 
the width, the parameter c and hence the corrections to the instanton density 

21 We used the criterium of assuming one standard deviation as x 2 = Xmin + b which is biased 
by the strong correlations between data for different momenta as just discussed. 
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estimates from the /c 4 -law in tab. |2J we will obtain roughly the same results. This 
could be expected after checking the rather similar behaviours of the function / 
obtained from both profiles in plot El This means consequently that, once some 
small width distribution is assumed and practically not regarding to the profile, 
the estimated instanton density is n ~ 6 — 8 fm -4 . 

All these estimates are anyhow in the right ballpark as far as different argu- 
ments [HI] seem to point towards an instanton density of a few fm -4 's. 

Discussing about realistic distributions. 

After estimating the effects of a small width distribution, we will discuss the 
consequences of some instanton distributions proposed in the literature. 

It is well known that the one-loop tunneling amplitude of classical minima 
in gauge theories gives the standard growth p(p) ~ p e for instanton radii distri- 
bution j2j. The classical interaction of instantons in the background is thought 
to introduce small-distance repulsion that strongly suppress large instantons and 
leads to a well-behaved partition function [Sj. Applying the Feynman variational 
principle to a liquid of BPST instantons, Diakonov & Petrov found indeed that 
this infrared growth is balanced as follows 



where p 2 = J °° dpp 2 p(p) and 7(72) = J °° dte t t n 1 is the standar Euler's gamma 
function and 



A very important remark is that eqs. (JBJ), particularized for BPST profile 22 and 
using j3{p) ~ P(p), combined with eq. leads to avoid the dependence on the 
density and to derive a prediction for the previously fitted parameter a-^pp, 



where we use j3(p) ~ 15 as in [3]. The agreement of such a prediction with our 
lattice estimate is remarkably encouraging since, for our non-BPST profile, this 
relation has been numerically checked and remain approximately valid. 

The goal of this section is anyhow to analyze the impact of the corrections 
to the equal radii approximation that we use to describe the low momentum 

22 In the BPST limit eq. © becomes eq. (QUJ where / dx 2 4> 2 = p 2 . 




(35) 




(36) 




(37) 
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behaviour of our lattice Green functions. For this purpose, we remark that the 
distribution given by eq. (|35|) is rather asymmetric around p, which is in practice 



usually defined as 



Its width, measured as 



p' 



2 

7 ' 



(3* 



leads us to conclude that the previous small approximation could have missed 
some not negligeable correction. This is why we compute in appendix [B] the 
Green functions and the MOM QCD coupling constant on the background of an 
instanton liquid with radii distribution given by eq. (|35|) and fit the latter to the 
power formula given by eq. (|31|). We do it for BPST and DP- inspired profiles 
and collect the results in tab. |3J As can be learned from the table, the main 
conclusions of the small width analysis remain unchanged except for the fact 



that the estimate of the density increases up to 
of the values of the parameter c. 



7 — 10 fm as a consequence 



Profile 


c 


e 


DPm 


1.55 


0.17 


BPST 


1.52 


0.25 



Table 4: Values of the parameters c and e in eq. Q31|) obtained in the appendix iBl for 
our DP-inspired and BPST profiles. 



Furthermore, the instanton radii distribution has to be varied altogehter with 
the instanton profile and density in the Diakonov & Petrov formalism. These 
authors found jlj: 

oil S(p) 

p(p) oc ^^exp(-(3n<T(p)p 2 >T(p)p 2 ) , (39) 

for the radii distribution of what they call "fremons" , i. e. the pseudo-instantons 
obeying the profile equation in (jZJ), where r(p)p 2 = 70/ dx 2 (p 2 (x/ p). s(p) is a 
function on p introducing a very small collective correction to the standard power 

P 6 - 

Two aspects which have lead Verbaarschot j2J an d Shuryak 122] to criticise 
the DP formalism are visible in the "fremons" radii distribution in eq. ()39|) : (i) 
the distribution concerns the radius, p, of one particular instanton immersed in 
some "bath" of instantons with average radius, p; how does one combine the two 
scales in the game, p and p, to define the instanton scale for (5 ? Why, for instance, 
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use (3(p) instead of P(i'pp) 1 ^ 2 ) ? (ii) Such an instanton scale is furthermore of the 
order of 0.3 fm and hence the coupling g in (3 is to be computed out of the 
perturbative regime! 

Any detailed evaluation of the impact of the "fremon" distribution implies 
answering these questions, in particular knowing the behaviour on p of (3{p) and 
how it modifies the gaussian in eq. (|35j). However, this is out of the scope of 
the present work. We adopt here a ^phenomenologicaV point of view: the DP 
formalism offers a large-distance damping for the profile that properly corrects 
the low momentum behaviour of Green functions and leads to a satisfactory 
description of our lattice data. The large- distance behaviour of the profile is 
governed by the parameter «dpP which is rather well estimated when we put 
P(p) ~ 15 in eq. (J37|) according to [3]. Moreover, had we assumed j3(p) ~ f3(p) 
and T(p) ~ 1 (this is strictely true for the BPST profile) we would recover the 
distribution in eq. (J35|) . 

As a final remark, if we try to compute the instanton density through any of 
the two eqs. (jlOI36j) by employing our best-fit parameters in the previous sec- 
tion, we will obtain ~ 2 fm -4 . Is something still missing in the DP approach? 
As a matter of fact, the third main aspect criticised by the authors above men- 
tioned E2J |^ is the too strong instanton interaction strength resulting from 
the DP formalism and it should be noted that a reduction of ^ by a factor three 
or four would allow all the pieces of the puzzle to match. Thus, our analysis 
detects somehow this effect pointed out in the above references regarding the 
instanton-antiinstanton interaction. It remains to investigate whether our results 
are compatible with those, for instance in [2*T] . where the Stream-line approach 
is used (the author found that 7q is about one order of magnitude smaller than 
the one from DP). 

We should recall however our "phenomenological" point of view: does the 
sum-ansatz of Diakonov & Petrov, after profile variation, take into account the 
effects included in the Stream-line analysis, except for some effective rescaling of 
7o ? A preliminar analysis of this question is avanced in appendix O where we 
show how one can build the gauge field through the Shuryak's ratio-ansatz |22j . 
that gives an approximative solution to the Stream-line equation, and match the 
sum-ansatz with our DP-inspired instanton profile in both near the center of any 
of the instantons and far from all them. This could explain why the independent- 
pseudoparticle approach given by the sum-ansatz accounts for the low momentum 
behaviour of the MOM QCD coupling constant and Green functions. This point, 
of course, deserves further attention. 

In summary, we dedicated this section to analyze the effects of radius dispersion 
on the /c 4 power-law that the instanton approach predicts for the MOM QCD 
coupling constant. Within the small width approximation we establish that the 
two main effects of radius dispersion appears in our fitting window as a reduction 
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of the power of k and as an overall pre-factor larger than 1. In particular, the pre- 
factor does not depend on the particular profile function and corrects the estimate 
of the instanton density from the MOM QCD coupling constant that becomes 
then compatible with those directly obtained from Green functions in the previous 
section, ~ 8 fm -4 . In fact, this instanton-density estimate is rather general and, 
as far as the low momentum behaviour of Green functions is dominated by the 
gauge field at large distances, is not only independent of the profile function 
but neither is affected by the sum-ansatz approach we used. Nothing seems to 
change by the use of realistic distributions. Futhermore, the analysis of radius 
dispersion in terms of realistic distributions gives us the possibility of connecting 

the "dynamical" parameter a^p-yp^ with only the instanton coupling throug (3 
and leads to an estimate of ~ 0.6 for the former, in remarkable agreement with 
our fits. 

6 Discussions and Conclusions 

In a previous paper |Hj, we had shown a k A dependence of the coupling constant 
in MOM scheme, which gave us a rather clean indication of the dominance of 
semiclassical solutions over low energy QCD dynamics, and some indication that 
these solutions might be instantons. 

In order to go further in the understanding of this very remarkable feature, 
we tried here a direct description of gluon Green functions in this framework. In 
order to achieve this, contrarily to what happens for the MOM coupling constant, 
the knowledge of the instanton profile is mandatory. The single-instanton BPST 
profile clearly fails to describe the low momentum behavior since it predicts a di- 
vergent gluon propagator for k — > in full contradiction with the lattice results. 
This is understood as an effect of the deformation of instantons far from their 
centers due to the influence of other instantons. Indeed, the BPST singularity of 
the gluon propagator is corrected when effects of instanton interactions are taken 
into account through a parameterization of the profile derived from variational 
methods [I], and it leads to a successful description of the low momentum be- 
havior of lattice two- and three-point Green functions (see fig. EJ). At this point, 
we cannot exclude other "instanton-like structures" as, for example, a significant 
amount of merons EI] which would of course modify our density estimates. 

We fit our lattice data for the Green functions in an instanton liquid model 
(ILM), below energies of ~ 0.9 GeV, and above ~ 0.4 GeV. Below ~ 0.4 GeV we 
have too few lattice data and furthermore the lattice computed gluon propagator 
gives a non null value for the k — > limit, in contradiction with the expectation 
of an ILM: 




(40) 
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Futhermore, the quantum corrections to the semiclassical solutions should be 
more visible when the latter vanish 23 . In fact, a signal that instanton approach 
could be missing some mechanism in the very low momentum regime is the fact 
that a very recent analysis within the SD approach leads to 24 

Gsd ~ {k 2 ) 2K ~ l , (41) 

with k ~ 0.595 [21] ■ An also recent SD-motivated work, based on rather general 
arguments, points towards k = 1/2 Last but not least, this region is expected 
to be strongly finite- volume sensitive and deserves a special study 

We have discussed the effect of the instanton radius dispersion. We find a 
subdominant but non negligible influence of the radius dispersion on the coupling 
constant: We find that a coherent description of both two- and three-point Green 
functions only emerges when the instanton radius reaches the vicinity of 1/3 fm, 
that is the value derived from phenomenological arguments. Furthermore, after 
invoking a radius distribution obtained from DP's variational met hods |4J, we 
correct the previously predicted fc 4 -power followed by a s {k) by a slight reduction 
of the power (~ k 3m ) and a 40 % increase of the prefactor. The fits of all our 
lattice data to a free-power law on k always produces a best-fit exponent lower 
than 4 but statistically compatible with 4 (We obtain 3.91(45) from a global fit 
combining all our lattice data sets). This seems to confirm the trend predicted 
from radius distribution, but the small number of points in the region of interest 
and the big errors lead to a large uncertainty that prevents us from a more 
conclusive statement. Once this radius dispersion is considered we get a good 
agreement of the density derived directly from a fit of the two- and three-point 
Green functions and from a fit to a s (k). 

The lattice spacings used in this study have been chosen rather large since we 
wished to reach low momenta with a not too large number of lattice points. It 
might be feared that these values are too far from the continuum limit, too much 
in the strong coupling regime, to be reliable. We did not see any sign of a non- 
smooth dependence of any quantity on the lattice spacing. However, for safety, 
it is advisable (and under progress) to follow-on these studies at, say, (3 = 6.0 
with as large a physical volume as used here. Larger statistics and larger physical 
volumes at a given lattice spacing are also needed. 

The lattice gluon propagator at k — > 0, has also open problems, for exam- 
ple in relation to the Zwanziger problem and critical exponent in relation 
to instantons. Other theoretical questions are pending such as the reason why 

23 In ref. the authors discuss about the impact of quantum effects on the large distance 
regime of instanton solutions. 

24 It is possible to find in the literature previous estimates for the critical cxponencnt, n, 
under different approximations' schemas as, for instance, k ~ 0.92 [2011231, « — 0.77 [58] or 

k= i nm. 
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quantum fluctuations only play a sub-leading role in this momentum range 25 . 

To briefly conclude, ILM can describe the low-momentum behaviour of gluon 
Green functions after modelling the instanton profile within the DP approach. 
Three are the parameters playing the game: a, the instanton radius and den- 
sity. The first one can be satisfactorily computed within the DP framework, the 
phenomelogical estimate of instanton radius leads to the best fits to lattice data 
and only the fitted density is much higher than its phenomenological estimate 
(although in agreement with another lattice estimates). Such a high instanton 
density leads to a large instanton packing fraction that makes hard a simple semi- 
classical approach to the gluon dynamics. However, instanton density should be 
reduced by including light dynamical quarks and it can not be then excluded that 
the semiclassical mechanism we use here might give account of the phenomeno- 
logical value. Our large packing fraction do not rule out a semiclassical approach 
to the full-QCD dynamics. 

Still we believe that we have given a series of rather convincing evidences 
of the influence of instanton-like structures on the low-energy QCD and more 
precisely on the low-momentum behavior of gluon Green functions. 

A Small and large momentum regimes of / 

The goal of this appendix is to perform the analysis, based on rather general 
grounds, of the small and large momentum regimes of the function / defined 
above in eq. ([29)1 . We will assume both regimes to be controlled by the Bessel 
function, J 2 (kpz), in the profile integral in that equation. 

For the small momentum case, the Taylor's series of the Bessel functions 
implies the following analytical expansion on kp, 

POO 

/ zdzJ 2 {kpz)(f){z) oc (kp) 2 (1 + \ p (kp) 2 + ...) , (42) 
Jo 

and we will then obtain 

«(*) = ^ (l + 6 ~f (48 + A8X p kY + 0(kY)) + O(^)) , (43) 

where, up to that order, only the coefficient X p keeps the profile function infor- 
mation. That coefficient is furthermore a negative one because of the alternance 
of signs in the Taylor expansion of J 2 , 

If we naively add quantum and classical contributions to the gluon propagator (Eq. 1)16(11. 
the perturbative quantum contribution would diverge. 
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and the correction to the /c 4 -power law is hence negative. However, this result 
is only general provided the analyticity for the kp expansion in eq. (}4*2*|) and, 
for instance, the BPST profile disregards this condition as can be seen by just 
expanding eq. (|3L)j) on kp, 



mm = {-^+ 2r ^ ln (f)} k2 p 2 + •■■ > ( 45 ) 

which is clearly non analytical in kp = 0. Nevertheless, we show in plot |H1 that 
f(kp) is negative for all k~p and the correction to the fc 4 -power law in eq. (}2*Bj) 
obeys to the same pattern as the one in eq. (f4l?)) with negative \ p . The reason 
for this behaviour of / can be mainly found in the large momentum limit. 

For the large momentum case, the oscillating nature of J2 leads the asymptot- 
ical behaviour of the profile integral in eq. (|4*2*|) to be dominated by the condition 
0(0) = 1 and hence by 



lim / zdzJ 2 (kpz) e~ az = -— . (46) 



Thus, for kp » 1, 



Pf(kp) 



- d 1 

V 



(kpf 



-2 . 



(47) 



As can be seen in plot El pf for both BPST and DP-inspired profile joins the 
horizontal asymptote —2 as soon as k~p ~ 7. 

Of course, this is not a rigorous proof of that corrections to the /c 4 -power 
law should be negative for any instanton profile (our fitting window is around 
k~p ~ 1), but leads us to reasonably assume that this is the case for most of the 
well-behaved profiles. 



B Some results for DP distribution of BPST in- 
stantons 

This appendix is devoted to present particular results for Green functions and 
for the MOM QCD coupling constant by employing the distribution for BPST 
instantons given in eq. (J35)) . We will give results for both BPST and DP-inspired 
profile functions. 
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The BPST profile. 

The use of BPST solution, <f)(\x\/p) = p 2 / (x 2 + p 2 ), implies neglecting the effect 
of neighbouring instanton's classical interaction. The IR behavior of gluon Green 
functions for this profile is: 



G^(k 2 ) 



k 2 ~ m n ( (5 



m 2 2m ~ 2 
Ak 2 



m/2 



n 



m 



6. 

/547T 2 

6"^ 



< p im I(kp) m > 



(48) 



'i + o(* 



x < 



2p 2 k 2 



for kp >> 1 



for kp <C 1 



where (3 = 6/g 2 at t * s the lattice parameter for the bare coupling and 



, .,„, mil + 2m) ( 3 , /9 + 2m\ , 2— 19 , \ 

U^ 2 ) = 1 [2jE ~ g + ^ ( v ^— J + ln 7^ - ln4 J • ( 49 ) 

Where 7^ the Euler constant ( 0.577216...) and ip(z) = r y'(z)/j(z) the Euler's 
"digamma" function. The first correction to the low momentum k A behavior of 
the quantity ((TJ), is then given by: 

121 , / 6-^,/ , Itk 2 22567\ \ , , 

valid for ~i?k 2 « 7/2. 

If we fit numerically the prediction of a s (k) to the power formula in eq. (|31|) 

(with p = \J~~p 2 = 1.5GeV _1 = 0.3 fm from table HJ) in the window (0.3-0.9) 
GeV 26 we will obtain: 



°-<*> = i^F (fc?) " 5 ' (51) 



26 This window is slightly larger than our lattice window (0.44-0.89) GeV. The question here 
is the validity of approximating cqs I|48I49|) by a power law which is shown to extend further 
than our lattice fitting window. 
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The DP-inspired profile. 

The same averages with our optimal profile parameterization in eq. (|13|) of 
section |U lead to the following low-momentum (k << odp) prediction for Green 
functions: 




G im) (k 2 ) - n [ J?- ^ VT 7( 5 + m) 



m 



(52) 



The result equivalent to eq. (|51|) for this profile, using the same fitting pro- 
cedure over the same window, with a phenomenological radius "p = 1.5 GeV and 
«dpP = 0.675 from table [TJ is: 

«•<*> = iS* • (53) 

C The ratio-ansatz 

In ref. j22] the following trial function, named ratio-ansatz, 



9 T? aa n a n 2 li^l 



*W = - , (54) 

i=I,A Ul 

where jji = \x — z % \, was proposed to avoid singularities not physically justified 
at the center of each instanton. In eq. 77° should be replaced by 77" when 
suming for anti-instantons as i = A. f(x) is a shape function that obeys /(0) = 1 
in order not to spoil the field topology at the instanton centres and that provides 
sufficient cut-off at large distances for the sum convergence. One of the physical 
motivations of such a large distances behaviour suggested in |22j is, in fact, the 
work of DP pE]. The author of [22], for simplicity, uses the following gaussian 
shape: 

f(x) = exp (-C^\ , (55) 

and obtain the coefficient C by the minimization of the action per particle for 
some statistical ensemble of instanton. However, the minimization of the repul- 
sion of instantons in matter leads, according to j3j , to the large distance behaviour 
given in eq. (|12|). Then, why not to use: 
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for x 







(56) 



Thus, if yi » pi for all % gB a behaves as expected following DP, 



9 b;( x ) ~ 2 £ 



i=I,A 




K 2 (anpyi) , 



(57) 



and as Hi/pi —> for any i = j 




- + K 2 (a DP y j ) 



K 2 (a D pyj) 



(58) 



In both large and small distances regimes we obtain the same through the sum- 
ansatz eq. (J2J) with the instanton profile eq. (|13|). The role of the coefficient C in 
eq. ()55|) is played by «dpP, and the minimization of the action per particle in [22] 
leading to the determination of C is acomplished by eq. (}3Tj) resulting from the 
DP variational procedure. Of course, we do not claim that the sum-ansatz with 
the profile eq. (J13)) leads to exactly the same results as the ratio-ansatz approach 
because the latter would imply that distances other than low and large ones 
have no influence on the Green functions. However, since the low momentum 
Fourier transform of the gauge field is dominated by large distances, at least the 
low-momentum behaviour of Green functions can be effectively described by the 
independent-pseudoparticle sum-ansatz approach. 
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